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Abstract 

Aghannan and Rouchon proposed a new design method of asymp- 
totic observers for a class of nonlinear mechanical systems: Lagrangian 
systems with configuration (position) measm'ements. The observer is 
based on the Riemannian structure of the configuration manifold en- 
dowed with the kinetic energy metric and is intrinsic. They proved 
local convergence. When the system is conservative, we propose a glob- 
ally convergent intrinsic reduced-observer based on the Jacobi metric. 
For non-conservative systems the observer can be used as a com- 
plement to the one of Aghannan and Rouchon. More generally the 
reduced-observer provides velocity estimation for geodesic flow with 
position measurements. Thus it can be (formally) used as a fluid flow 
soft sensor in the case of a perfect incompressible fluid. When the cur- 
vature is negative in all planes the geodesic flow is sensitive to initial 
conditions. Surprisingly this instability yields faster convergence. 

Keywords Riemannian curvature, geodesic flow, non-linear asymptotic 
observer, Lagrangian mechanical systems, intrinsic equations, contraction, 
infinite dimensional Lie group, incompressible fluid. 

There is no general method to design asymptotic observers for observable 
non-linear systems, but some specific types of linearities have been tackled 
in the literature. In particular work has been devoted over the last few 
years to observer design for systems possessing symmetries. Notably the 
case of a finite-dimensional group of symmetries acting on the state space 
[21 El [131 E] , and left-invariant dynamics on a Lie group have been considered 
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[H] . Symmetries correspond in general to invariance to some changes of units 
and frame. The invariance to any change of coordinates was raised by [3] who 
designed an intrinsic observer for a class of non-linear systems: Lagrangian 
systems with position (configuration) measurements. The aim is to estimate 
the velocity, independently from any nontrivial choice of coordinates, and of 
course never differentiate the (noisy) output. The observer was adapted to 
the specific case of a left-invariant system on a Lie group by [T^. Observer [3] 
is based on the Riemannian structure of the configuration manifold endowed 
with the kinetic energy metric. This geometry had already been used in 
control theory of mechanical systems (see e.g. [HI [ID])- The convergence of 
the observer is local. 

According to the Maupertuis principle, the motion of a conservative La- 
grangian system is a geodesic flow for the Jacobi metric, intrinsically defined 
after the kinetic and potential energies, up to a time reparametrization. In 
this paper we consider the general problem of building a reduced observer 
to estimate the velocity for geodesic flow on a Riemannian manifold (mo- 
tion along a geodesic with constant speed). Under some basic assumptions 
relative to the injectivity radius (also formulated in [3]) we have the follow- 
ing results (theorem 1). If there is an upper bound A > on the sectional 
curvature in all planes, the reduced observer is exponentially convergent (the 
convergence is global in a sense to define) as long as the gain is bounded from 
below by a linear function of \/A. Unfortunately the higher the gain is the 
most sensitive to noise the observer is. An even better situation occurs when 
the sectional curvature is non-positive in all planes: the reduced observer is 
globally exponentially convergent, without restrictions on the gain. In fact, 
for a fixed gain, the more negative the curvature is the faster the observer 
converges. This feature is surprising enough as negative curvature implies 
exponential divergence between two nearby geodesies, and thus "amplifies" 
initial errors. This is a major difference with [3] who used additional terms 
precisely to cancel the effects of (negative) curvature. 

For mechanical Lagrangian systems the observer of Aghannan and Rou- 
chon is only locally convergent. In the absence of external forces the reduced 
observer offers a good alternative, as it is globally convergent. When there 
are external forces, the reduced-observer can be used as a complement to [3] . 
The gain must be chosen large enough so that the reduced observer converges 
before the energy varies significantly. If so, it provides the observer [3J with 
an initial estimated velocity close to the true one. 

The reduced observer is also applied, formally, to a basic velocimetry 
problem: compute the velocity of a perfect incompressible fluid observing 
the fluid particles. The principle of least action shows that the motion of an 
incompressible fluid can be viewed as a geometric flow. We consider the case 
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of a two-dimensional fluid. As the convergence properties of the observer 
depend on the sign of the curvature, we will use results and heuristics of 
Arnol'd [5]. Following them we show theoretically that global convergence 
could be expected for a large class of trajectories, since the curvature is 
positive only in a few sections. This latter fact also implies instability of the 
flow, and Arnol'd gives an interpretation relative to the difficulty of weather's 
prediction. Note that the problem tackled is nontrivial, as the system is 
nonlinear, infinite dimensionnal, and possibly sensitive to initial conditions. 

In Section 1 we give the general motivations introducing Lagrangian sys- 
tems on manifolds and Maupertuis' principle. In Section 2 we introduce the 
observer. In Section 3 we consider applications to some mechanical and hy- 
drodynamical systems. In Section 4 the convergence in the case of positive 
constant curvature is illustrated by simulations on the sphere. 



1 Lagrangian systems on manifolds 

Consider the classical mechanical system with n degrees of freedom described 
by the Lagrangian 

= -9ij{q)q'q^ -U{q) 

where the generalized positions q & M are written in the local coordinates 
(?')j=i...n, g{,q) = (ygij{q))i=i...n,j=i...n is a Riemaniann metric on the configura- 
tion space M, and ?7 : i— >■ M is the potential energy. The Euler-Lagrange 
equations write in the local coordinates 

One can prove [1] using ^§^gjk = ~g^^^§^ where are components of g^^ 
that ([T]) writes 

q' =-^M)i'i' + wU (2) 
A curve 7(t) which is a critical point of the action 



where the Christoffel symbols F*^ are given by F*^ = + -q^ 



s{^)= r Li^{t),^{t)dt) 

Jo 

among all curves with fixed endpoints satisfies the Euler-Lagrange equations 
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1.1 Lagrangian system in a potential field 



Consider a conservative Lagrangian system evolving a an admissible region 
{q & A4 : U{q) < E}. The energy of the system E = T{q,q) + U{q) = 
\gij{q)q^q^ + U{q) is fixed. According to the Maupertuis principle of least 
action (see e.g. p]), in the Riemannian geometry defined by the Jacobi 
metric g'^^{q) = 2{E — U{q))g^\q) and the natural parameter r such that 
^ = 2{E — U{q{t))), the geodesic fiow is a solution of the equation of motion 

([2]). Indeed if the f*;, are the Christoffel symbols associated to the metric g 
we have 

|i«- + f;.w|:«'|:«' = o (3) 

which writes intrinsically 'V dq_jr = and defines the geodesic fiow (V is the 



dT 



Levi-Civita covariant differentiation of the Jacobi metric). 



1.2 Geodesic flow and holonomic constraints 

A material particle constrained to lie on a manifold moves along a geodesic 
[5]. Indeed E = T, g = 2Eg, dr = 2Edt ensure the energy T is fixed. Ac- 
cording to Maupertuis' principle the motion minimizes \Jgijjpq^jpq^dT = 

{l/y/2E) \J d^ili^^li^^^^ which is proportional to the geodesic length in 
the metric g. More generally an inertial motion of a Lagrangian system with 
k holonomic constraints can be viewed as the inertial motion of a particle 
constrained to lie on a submanifold of dimension n — k (see e.g. |5] p 90). 
A conservative Lagrangian system in a potential field with holonomic con- 
straints satisfies the Maupertuis' principle on the configuration submanifold 
of dimension n — k. 



2 An intrinsic reduced-observer 

Let us build an observer to estimate the velocity g of a point q moving along 
the geodesies of ^A with constant speed, when the position q is measured 
(with noise). First suppose M. = endowed with Euclidian metric. Let 
q = V and v = 0. Consider the following reduced-observer 



It can be interpreted as a kind of "elastic" trailer hooked up to q. It allows 
to estimate the position and speed of the system 

^ = q — \v (and thus ^ = v) (5) 

We want to prove i — ^ ^ and -0 — f where v = We have 

V = —{1/ \){v — v). If A > 0, -0 ^ f since f is a constant vector. But 
C, = q — Xv, thus — ^ and asymptotically ^ is moving behind q at 
distance A||f ||. Now if Ai is any Riemannian manifold consider 

d 



j^i = -^grad^D\i,q) A>0 (6) 

where D{C,,q) is the geodesic distance between ^ and q. The dynamic does 
not depend on any choice of local coordinates in M", and is a generalization 
of (jl]). We want to prove that D{^,^) — > where ^ is a point following q at 
distance A||w|| on the geodesic {q(t) : t > 0}. The parallel transport 

of to the tangent space at q along the geodesic joining ^ and q is an 
estimation of f = g. With the following definition, v is no more noisy than q 

Theorem 1. Let A4 be a Riemannian metric. LetT < oo. Lett i-^ q{t) G Ai 
satisfy V^g = for t G [0,T]. Let ^(t) = expg(^^^{—X\\q\\). Consider the 
observer Let ([H]) be the condition 

Dm,m < e--^'Dim,m) vt g [o,r] (s) 

• Suppose the Riemannian curvature is non-positive in all planes. If for 
all t G [0,T], D{^(t),q{t)) is bounded by the injectivity radius I{t) at 
q{t) (i.e. there exists a unique geodesic joining ^ andqit)), ^ is true 
for all A > 0. When the manifold is complete and simply- connected 
(Hadamard manifold), the injectivity radius is infinite ( Cartan-Hadamard 
theorem) and ([H]) is always true. In particular ^{0) can be chosen ar- 
bitrarily. Moreover 

x\m - m\\ < Dm,m) < e--^'D{m,m) vt g [o,t] (9) 

• Suppose the sectional curvature is bounded from above by A> 0. Take 
X > 4||^|^^ - & is true as long as the distance D{^{t),q(t)) remains 
bounded by max(^^,/(t)) for all t G [0,T]. // the manifold is simply 
connected, ((H]) is true as soon as i5(^(0), g(0)) < 
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Under the assumptions of theorem 1 (second case), if the observation 
noise amphtude is less than the observer is also globally convergent in 

the following sense: take ^(0) equal to the initial position measurement. We 
have (IE]) without any information on the initial velocity. 



proof The proof consists in two differential geometry lemmas. 

Lemma 1. Let M. he a smooth Riemannian manifold. Let P E M. he fixed. 
On the suhspace of M. defined hy the injectivity radius at P we consider 

j^x = -j^^d., D\P,x) A>0 (10) 

// the sectional curvature is non-positive in all planes, the dynamics is a 
contraction in the sense of IT^ . i.e, if 6x is a virtual displacement at fixed t 
we have 

^^\\8x\\l < ~\\6x\\l (11) 
where \\ \\g is the norm associated to the metric g. 

If the sectional curvature in all planes is hounded from ahove hy A > 0, (ITT]) 
holds for D{P,x) < ^l/{A^/A). 



proof The virtual displacement is defined [12] as a linear tangent differen- 
tial form, and can be viewed by duality as a vector of TMI^,. Let us define a 
surface S. Let 70 be the geodesic joining P to x. A point which is infinitely 
close to X in the direction Sx is linked to P by a geodesic 70 + ^7. The 
directions defined by this latter geodesic and 7 at P span a 2-plane tangent 
at P. All the geodesies having a direction tangent to this 2-plane at P span 
a smooth surface S embedded in which inherits the Riemannian metric 
g. S is invariant under the flow ffTUl) . as the gradient term is tangent to the 
geodesies heading towards P. We have x G S and 6x G T^S C T^A^. 

Following [18j (p 177) we use specific coordinates on S called "polar 
coordinates". Let ei, 62 be an euclidian frame of TpS for the inherited metric 
and ei be tangent to 70. We define $ : {<J,0) expp{a cos 6ei + cr sin 6^62). 
S is parameterized by a, the geodesic length to P, and 9, the angle in TpS 
with ci. In the polar coordinates, the elementary length is given by 

ds'^ = da'^ + G{a, e)de'^ 

and G satisfies the initial conditions \/G = and = 1 at cr = 0. Ac- 
cording to a classical result the Gauss curvature at the point (cr, 6) is given 



6 



by K{a,9) = ^ ^^2'^'^'' ■ We will prove (lemma Ej) that the Gaus- 

sian curvature at m = $((T, 6') G S is bounded from above by the sectional 
curvature in the tangent plane to S at u: K{a,9) < KseciTyY,). Suppose 
KseciTyY,) < 0. It implies K{a, 6) < 0. Along 7 we have 



da"^ da^ ^ ' ' da ' \^ da 



(12) 



and thus ^(c"|^) > |^ which yields by integration a^ > G since G{a^ 0) = 
0. In the polar coordinates the dynamics (|TU|) reads 

a = -\a] e = {] 
A 

Writing ||5x||2 = a^da"^ + l3'^G{a,Q)56'^ we have along the geodesic 70 (pa- 
rameterized by a and ^ = 0) the following inequality, proving (fTTl) . 

Usxf = -2^5a^ - 2P^^^^p^59^ < -hsxf (13) 
dt A A ua A 

Suppose now KseciTuT,) < A. Let us find an upper bound on a under 
which ( ITTl) holds. Let x{a) = a/ G{a, 0). It suffices to find an upper bound 
above which (1121) is not true, or a such that x'{a) = ^/Ax{a). We have x" = 
—Kx, x(0) = 0, a;'(l) = 1. Let y{a) = sm.{\/~Aa) / \/~A. A Taylor expansion 
in shows that there exists a > such that x{a)/x\a) < y{a)/y'{a) for < 
a < a. Let (Pk{<^) be the polar angle of the point {x'{a),x{a)) in the phase 
plane, and ^A{cr) the polar angle of (cos(a/]4(t), sin{\/~A.a) / ^/A) . It is shown 
in the proof of the Sturm Comparison theorem in [3] that ^Pa{o') > ^k{(^) 
(as long as cr < ti/\/~A). Thus as long as \/~Aa < arctan 1 i.e. a 
we have x'{a) > ^/Ax{a). Note that we proved at the same time Tr/(A^/A) 
is a distance under which the exponential is an injection. If the manifold is 
simply connected, the exponential is a surjection and tc/(A\^) is less than 
the injectivity radius. 

Lemma 2. Let M. he a smooth manifold. Let P E M.. Let E he a two- 
dimensional vectorial space of TpM. . Let uj he a neighhorhood of P in E 
such that the restriction of the exponential map a to u is a diffeomorphism 
in M.. S = aiuo) is suhmanifold of dimension 2. Its Gaussian curvature at 
any u & is hounded from ahove hy the sectional curvature in the tangent 
plane to S at u: 
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proof The proof is written after computations of Ivan Kupka. Let S be 
the surface constructed at lemma [T] and $ the map associated to the polar 
coordinates. The metric inherited on S writes 

= 11^^^ + ^^^"^ = "l^ll^^^ + 2 < ^' ^ >^ ^^^^ + "l^ll^^^ 

For fixed 6 the curve 70 : cr i— ^(a,6) is a geodesic and thus |||f ||g = 1- 
Let J{o-,6) = ||. For fixed 9, Jg : a ^-^ J{(^-,9) is a Jacobi field along 7e. 
Moreover J(0,6') = — sin^ei + cos6'e2 is orthogonal to this geodesic in P. 
It is well known (Jacobi field properties) that it implies J is orthogonal to 
761 at any point. Thus < |f , || >g= and ds^ = da"^ + \\J{a,6)\\'^gd6'^ , and 
the Gaussian curvature is given by K{a,9) = ||j(~g)|| ^ ^^-^^^f^^^" (see [IE])- 
Consider the Levi-Civita covariant differentiation V of the metric g. We 
have ^^Mk = <^-f-'^)J-'^)> and 

da \\J{a,9)\\g 

d^J{a,e)\\g _ <VlJ{a,e),J{a,e)> <V^J{a,e),V^J{a,e)> 
da' \\J{a,e)\\, + \\J{a,e)\\, 

_ < V^J{a,9),J{cr,9) >' 

According to the Jacobi equation we have J((t, 9)+R{J{a, 9), ^^|^)^^|^ 
0. Thus the Gaussian curvature of S satisfies 



Kia,9) = Ksec{Tu^) + 
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< V^Jjcr, 9), J{a, 9) >' -\\J{a, 9)\\l\\V,J{cr. ^^"^ 
\\Jia,9)^ 



where Ksed-l-u^) = \\j{ae)p sectional curvature at 

u = $(cr, 6'), and where we used that J(cr, 9) is orthogonal to Cauchy- 
Schwarz implies that the fraction above is negative and K{a, 9) < Ksec{Tu'^)- 
We can now proove theorem 1. Suppose the sectional curvature is nonpos- 
itive in all planes. Using the contraction [12] interpretation in the appendix 
of [3] we see that if ^1,^2 are solutions of ([HD which remain at a distance of 
q bounded by the injectivity radius for < t < T we have 

DiUt),Ut)) < e-i*D(a(0),6(0)) Vt G [0,r] 

since by lemma 1 ([6]) is a contraction. The system "forgets" its initial 
condition. So ([8]) holds if ^(t) is a solution of (Q. This is true since 
= iD{^,q) = \\v\\, - WiWg = \\v\\g - {Di^,q). To prove (0) notice 

^ implies G{a, 9)d9 > a'd9. Thus the angle a between 
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Figure 1: Left: Geodesic deviation on a manifold of negative curvature. (fTOj) 
writes in polar coordinates {cr,9): a = —j;6 = 0. The distance \\Sx\\ between 
neigbors Xi and X2 decreases at a rate at least j. Right: Geodesic deviation 
on the sphere. 

V and V in TgAi is smaller than the angle /3 corresponding to the Euclidian 
case K = 0, G{cr,6) = cr^. In this latter case q,^,^ is an Euclidian triangle 
and thus \\v — q\\ < D{^,^)/X (make a drawing to see the homothety). 

Suppose the sectional curvature is bounded from above by A and A > 
Then if Dim, ?(0)) <^^e have q{t)) < ^ for all t > 

0. Indeed iDii,q) << giad^Di^, q),v > -WfJUa < hh " jDiiq). Thus 
-D(^)Q') = implies ^D{^,q) < 0. By lemma 1 ([H]) is a contraction. ([H]) 
holds by the same token as last paragraph. When the curvature is constant, 
it is easy to prove that \\v — q\\ — >■ exponentially at rate 1/A since G{(J, 6) = 
^^^^^ ■ Intuitively we suspect it remains true if the curvature is not varying 
too fast along the geodesies but we have no proof for this claim. Nevertheless 
for fixed q, ^ and time t, if D{^, ^) < vr/ (Ay/A) we have in the limit H'W— ^|| — > 

asD(e,o--o. 



3 Applications 

3.1 Lagrangian mechanical system 

Proposition 1. Consider any Lagrangian system in a potential field in the 
admissible region defined by U < E. The observer 

jj = -\iE - U{q)) D%q) (14) 



is such that the theorem 1 is valid in the Maupertuis time and Jacobi metric. 

proof To prove this proposition, one can apply the Maupertuis' principle 
(see section [TTTll . In Maupertuis' time r = 2(E — U{q{t))dt the motion is a 
geodesic flow on the configuration space with modified metric g, with \\v\\g = 
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1 as ^ 2{E - U)g,,iq)^^i§y = 1. The observer defined 

by = —■^grad^Dg'^(^,q), A > where Dg is the distance associated to 

Jacobi metric, is such that v = 'T^/^^^^C is an estimation of ^g. 

For instance, R. Montgomery studied in a recent paper [T3] the Newto- 
nian equal-mass three bodies problem, with zero momentum and when the 
potential is taken equal to the Jacobi metric has negative curvature ev- 
erywhere (except at two points). The reduced observer (HM is thus globally 
convergent for a three bodies system which is sensible to initial conditions. 

Remark 1. For a conservative system, the total energy E, or equivalently 
the initial kinetic energy needs to he known to compute the Jacobi metric. 
However no information about the initial direction of the velocity is required. 
In the particular case of an inertial motion, the metric g is used and no 
information on the initial velocity needs to be known ( except an upper bound 
to tune the gains in the case of positive curvature). 




Figure 2: Ball and Beam 

Remark 2. Let us consider now a non- conservative system: the ball and 
beam of ^3] with a torque control u (see figlE)- The observer is only locally 
convergent. Observer ([6]) can be used complementarily to provide a globally 
convergent estimator with the following little experiment. Set ^(0) = q'(O). 
Maintain m = (no control). The characteristic time of convergence of the 
observer ^ is tq = \ in the Maupertuis time. After a few tq the observer 
([6]) provides the observer of Aghannan and Rouchon an initial estimation of 
the velocity close to the true one and from that moment u can vary freely 
again: [3] converges. There are two restrictions: the initial energy must be 
approximately known, and the speed must be bounded from below during the 
experiment. 

3.2 Motion of a perfect incompressible fluid 

The goal of this section is to show that the reduced observer could possibly 
be applied to more complicated systems. No formal proof is given but only 
heuristic discussions. The observer could be used in velocimetry as a (soft) 
velocimeter when the fiow is seeded with particles that can be observed and 
is modeled by Euler equations. 
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3.2.1 A reduced observer 



Let us first introduce some results and notations of [5], EJ |T7j. Let f2 be a 
domain of bounded by a surface SQ. Let v be the velocity field of an ideal 
incompressible perfect fiuid with density p which fills the domain Q. The 
motion is described by the Euler equation 

+ {v.\/)v = -^\/p (15) 

where p is the pressure. Let SDiff Q be the Lie group of all diffeomorphisms 
that preserve the Euclidian volume. Its Lie algebra U is the set or all vector 
fields of fl of null divergence, and tangent to 6fl. Consider the scalar product 
on the Lie algebra 



yv,wEU, <v,w>=p / / / v{x).w{x)dx (16) 




Let v{t) G W be a solution of ffT^ . Let ^^(a;) = z{t) be the position of every 
fiuid particle after a time t, obtained by integration on [0, t] of the system 
-^z = v{s,z), z{0) = X. (pl is a diffeomorphism for any t > 0, and the 
motion of the fiuid is described by a curve t i— 0J' on SDiff VL. Suppose t is 
fixed. After a small time r the diffeomorphism describing the fluid will be 
expj^(r'y(t))0j up to second order terms in r. It implies v{t) = DR^^vyi-j^cjf^ 
where DRg denotes the tangent map induced by right multiplication by g 
on the group. Thus the kinetic energy of the fluid T = | > deflnes 

a right-invariant metric. The least action principle implies that the fluid 
motion t i-^ 0^ is a geodesic flow on SDiff VL endowed with the kinetic energy 
metric. Thus V^^v = where the Levi-Civita covariant differentiation V^*" 
is given by V^^ff = + (t? ■ V) ?f + Va and a is a real function such that 
V^*" ff eU. For flxed t, deflne the virtual displacement 50(s, ■) corresponding 
to 6x in lemma 1: 

50(s,x)^0f^'^(x)-0^» 

|17j proves it can be identifled to an element of U and it satisfles the Jacobi 
equation (V^'")^50 + A^^Scj)) = where < v4^(rf), 17 > is the curvature of the 
plane section spanned the orthonormal frame v, ff. 

The reduced observer is deflned intrinsically and can formally be applied 
to this fluid velocity estimation problem. The theorem 1 is valid, as the proof 
is only made of intrinsic calculations, and its core is the Jacobi equation which 
gives conditions under which cr-^G > G. The observer's state ^ is a virtual 
fluid deflned as a solution of where q is replaced by (p^. Using the right 
invariance of the metric, one can consider ({t) = ^[t) o at any t > 0, 
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and the update is given hj v = —DR^vW G T^vSDiS Q where ( = (f)i. Thus 

only the inverse of the exponential map (p^ is required. Note that ( must 
remain in the identity connected component. 

3.2.2 Discussion on the convergence and curvature 

When the curvature is bounded from above by —B^ < 0, the geodesic flow is 
sensitive to initial conditions, and amits ergodic properties [5]. Surprinsingly, 
in this case the observer is globally exponentially convergent by theorem 1. 
When there are always sections with negative curvature along a geodesic, it 
is commonly assumed that the sensity to initial conditions is still valid. 

We have the following formal convergence result: consider a sinusoidal 
parrallel stationnary motion of a fluid in the tore = {{x,y), x mod 27r, y 
mod 27r} given by the current function = cos{kx + ly) with A;, / G N, and 
the velocity vector field v = rot ip. Take ^(0) = 0o &■ Then ||'0(t) — v\\ 
converges exponentially to 0. The proof is obvious as both points belong to 
the same geodesic. But one can expect a great robustness to measurement 
noise. Indeed [5] proves the motion defined by is a geodesic of SDiff T^, 
and the curvature is non positive in all planes containing v{x, y). Moreover it 
is zero only in a family of planes of null measure. But by theorem 1 negative 
curvature implies global stability, and small positive curvature implies a large 
basin of attraction. 

More generally, Arnol'd [5] considers the group SoDiff of diffeomor- 
phisms preserving the center of gravity. Calculations show the curvature is 
positive "only in a few sections". He suggests, as a concluding remark of 
appendix 2, to consider the mean curvature along paths to characterize the 
stability of the flow. As a consequence, if the atmosphere was a bidimensional 
incompressible fluid on the earth viewed as (identify opposite sides of the 
planisphere), the wind should be known up to 5 decimals for a two-months 
weather's prediction. Following this suggestion, as the curvature is positive 
in only in a few sections, one could expect a good global behavior of the 
observer. 

Remark 3. One could object that the problem is as simple as identifying 
a geodesic from position measurements. But trying to fit the geodesic by a 
least-square type method (for instance) would be less robust to perturbations: 
a slight ponctual change in the velocity angle at some point leads to large 
changes in the trajectory. However the observer follows (p^. Moreover the ob- 
server is intrinsic, its gain is easy to tune, and no differentiation is required. 
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Figure 3: Simulations on the sphere for A = 7r/4. Left: q (plain line) et ^ 
(dashed line). Right: velocity v (plain line) and estimation v (dashed line) 



4 Simulations on the sphere 

Consider the inertial motion of a material point constrained to lie on the 
sphere S^. The speed is constant (see section \UI\\ and assumed to be equal 
to 1. One can always choose coordinates q = {qx,qy,qz) G M.^ such that the 
motion writes: qr,{t) = cos(t), qy{t) = sin(t), q^{t) = 0. Let ^ = {C.x,^y,^z) e 
R^. The observer equation ([6]) writes 

II (gAOA^Il' 

±i = 1^ ((gApAq 

dr' y II (gAOA^ir 
±, ^ 1 ((gAQ AQ. 

fit a'^ II (gAO aIII 

where A < and ip is the angle between q and ^. As the geodesies of the 
sphere are great circles, ip is the geodesic length between those two points. 
The inital conditions are : g(0) = [1,0,0]^ and ^"(0) = ^[0,1,1]^. To 
simulate the sensor's imperfections a white noise whose amplitude is 20% of 
the maximal value of the signal was added. ^ converges to the equator, and 
asymptotically follows g at a distance |A|. The parallel transport v of 
is an estimation of f . In fact for A < 7r/2 the observer always converges in 
simulation, and for A > tt/2 it does not (see fig [I]). 
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5 Conclusion 



We designed a nonlinear globally convergent reduced observer for conserva- 
tive Lagrangian systems. The observer is intrinsic and converges despite the 
effects of curvature: instability of the flow and gyroscopic terms. The tuning 
of the gains is simple. The only gain is a scalar to be chosen in function of 
the noise and the maximal curvature. It can be used for velocity estimation 
for all systems described by geodesic flows (VvV = 0), notably conservative 
Lagrangian system, and the motion of an incompressible fluid. Using the 
Maupertuis principle this work could be extended to the case of a mixture 
of compressible fluids [16j . 

Unfortunately when the motion is described by V^v = S{q) with S known 
(Lagrangian system with external forces) the reduced observer does not con- 
verge. Including such terms S remains an open question. As a concluding 
remark, note that the article gives insight in the link between convergence 
and geometrical structure of the model in the theory of observers (comple- 
menting the work of [3] and more recent results [S]). The observer's principle 
can be viewed as the simplest pursuit algorithm on manifolds (proportional 
feedback). A quantitative link between the effectiveness of the pursuit and 
the topography of the environment (Riemannian curvature) has been estab- 
lished. 
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